ANALYSIS 55(282), June 1, 2019 


= DISCOVERY 


? 
, : 


Performance study of three-phase induction 
motor driving a load 


Oti Stephen Ejiofor’, Nwosu Cajethan Abuchi*, Nnadi Damian B N°, 
Ogbonnaya | Okoro* 


'Department of Electrical Engineering, University of Nigeria, Nsukka, Enugu State Nigeria. E-mail: stephen.oti@unn.edu.ng 
*Department of Electrical Engineering, University of Nigeria, Nsukka, Enugu State Nigeria. E-mail: cajethan.nwosu@unn.edu.ng 
3Department of Electrical Engineering, University of Nigeria, Nsukka, Enugu State Nigeria. E-mail: damian.nnadi@unn.edu.ng 
4Dept. of Elect./Electronics Engineering, Michael Okpara University of Agriculture, Umudike, Abia State, Nigeria. E-mail: 
oiokoro@yahoo.co.uk 


Article History 

Received: 11 April 2019 
Accepted: 22 May 2019 
Published: 1 June 2019 


Citation 
Oti Stephen Ejiofor, Nwosu Cajethan Abuchi, Nnadi Damian B N, Ogbonnaya | Okoro. Performance study of three-phase induction 
motor driving a load. Discovery, 2019, 55(282), 279-290 


Publication License 


(+) 
ey © The Author(s) 2019. Open Access. This article is licensed under a Creative Commons Attribution License 4.0 (CC BY 4.0). 


General Note 
Article is recommended to print as color digital version in recycled paper. 


ABSTRACT 


Induction motors are the most indispensable and widely used electric motors in industry. During start-up, loading and unloading 
operations, the asynchronous motor draws large currents, produce voltage dips and oscillatory torques. It is therefore important to 
model and study the asynchronous machine in order to predict these phenomena. Lots of work has been done on induction motor, 
ranging from transient and steady state responses to stability and performance studies. However, the effects of temperature 
variation during motor operations, variations in external rotor resistance and variable frequency drives have not received the desired 
attention. Hence, this work presents a step by step Matlab/Simulink implementation of a three-phase induction motor driving a 
load. Herein, a model of 126hp rated three-phase induction motor driving step input loads at different levels was demonstrated. 
Motor parameters were varied in order to study the effects of such variations on transient response and general performance of the 
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motor. The electromechanical torque and motor speed responses of the motor were examined. Consequently, recommendations for 


optimal motor operation were proposed. 


Keywords: Asynchronous, Matlab/Simulink, Transient, Performance, Load, Step input 


1. INTRODUCTION 


The three-phase induction motor is the most common of all industrial motors. The use of induction motors has increased 
tremendously since the day of its invention. They are being used as actuators in various industrial processes, robotics, house 
appliances and other similar applications. The growing popularity of the individual drive system in small workspaces, and the 
increasing use of portable electric tools, has stimulated development of motors, which, though of very small size, are of excellent 
electrical and mechanical design, and form miniature power units of great reliability and good performance [1]. 

During start-up, loading, unloading and other severe motoring operations, the asynchronous motor draws large currents, 
produce voltage dips, oscillatory torques and can even generate harmonics in the power system [2]. It is therefore important to be 
able to model the asynchronous machine in order to predict these phenomena. 

In the past, researchers have developed their own software packages for dynamic modeling of induction motor [3]. It is 
unnecessary to develop user-written software for dynamic model of induction motor when proprietary software package such as 
Matlab/Simulink, licensed by MathWorks which makes simulation design easy and more efficient is easily available. 

Though some researchers have done similar work to this topic [2, 4, 5,6, 7,8], most of them simulated the induction motor model 
under no load condition while others with a specific load. The induction motor static model as developed in [4] did not take into 
consideration the effect of the load inertia. The load inertia is a very significant factor when it comes to stability and motor 
performance studies [5]. While [6] did a work on induction generator operating at a speed of 6600rpm, [7] studied the induction 
machine driving a pump load considering the steady and transient states in the analysis. Computer simulation of induction machine 
coupled to a mechanical load was well analyzed in [8] but in [9], 7.5kW squirrel cage induction motor used for study was modelled 
using Simulink block and the work neglected number of poles in rotor speed equation, it varied only the stator resistance in the 
motor performance study. Presently, this work is aimed at modelling and simulation of 126-hp rated, three-phase induction motor 
driving different load levels of ONm, 100Nm, 550Nm and 1200Nm using Simulink blocks. Furthermore, the stator and rotor 
impedance of the motor were also varied so as to achieve optimal operation which guarantees an optimal use of the three phase 
induction motor in different industrial applications. It is known that the speed of rotation of the d-q axis can be arbitrary although 
there are three preferred speeds or reference frames as presented in [2]: 

(a) The stationary reference frame where the d-q axes do not rotate. 

(b) The rotor reference frame where the d-q axes rotate at rotor speed. 

(c) The synchronously rotating reference frame where the d-q axes rotate at synchronous speed. 


Meanwhile, the following assumptions as noted in [9] were made to derive the dynamic model: 
i Uniform air gap. 

ii Balanced rotor and stator windings, with sinusoidally distributed mmf. 

iii | Saturation, eddy current and temperature effects are neglected 

iv Skin-effect neglected 


Direct quadrature transformation of three phase induction motor 
In balanced three-phase circuits, application of the dq transform reduces the three AC quantities to two DC quantities [11]. The two 
stages used to achieve dq transformation are the Clarke's and the Park's transformation. 


Clarke's Transformation 
Three-phase voltages varying in time along the axes a, b, and c, can be algebraically transformed into two-phase voltages, varying in 
time along the axes a and B by the following transformation matrix [12]. 


Fa] _ 4 =172 =41724/7¢ 
ire =2/3|° 3/2 —V3/2 ty nadandanedtemetianaadhateaean lies (1) 


www.discoveryjournals.org OPEN ACCESS 


Page2 OU 


ARTICLE 


The inverse transformation can also be obtained to transform the quantities back from two-phase to three-phase: 


Fa 1 0 
Fb| =2/3|-1/2 V3/2 ial Para iene aeenacneeers (2) 


Fc —1/2 3/2 


Park's Transformation 
The two-axis orthogonal stationary reference frame quantities are transformed into rotating reference frame quantities using the 
park's transformation. This is expressed by the following transformation matrix: 


[ral 2 cosO | & 
Fq —sin® cos@l [FB | 


The inverse transformation can also be obtained to transform the quantities back from two-phase rotating to two-phase 
stationary frame as follows [13]: 


a 2 Gee ee EB % 
FB sind cos@ | [Fg | 


Dynamic Model equations for 3-phase induction motor 

The dynamic model of induction motor in arbitrary reference frame can be represented by using flux linkages as variables. The 
model equations can be generated from the dq equivalent circuit of the induction machine [11, 14]. Figures 1 and 2 below show the 
equivalent circuit of the induction machine. 


ids Rs Lls=Ls-Lm Llr=Lr-Lm Rr. ide 
ids . i: at 


Vds Vdr 


Figure 1 d-axis equivalent circuit of the Induction Motor 


igs Rs Lis =Ls-Lm Llr=Lr- Lm Rr. 4 


Figure 2 q-axis equivalent circuit of the Induction Motor 
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Flux Linkage Equations 


The flux linkage equations waa with this circuit can be found as follows [2]: 


“« = wb [vas — — Ads += a ~ (Amq — aqs)|..... ee one (5) 
ee = wb [vas — “Ags + — (Amd — Ads) |. ... ene © 
da (we-wr) R 
7 = wb [var — ae dr + mH (Amq — har) -dbs Seeeeiaieaee (7) 
Be = wb [Var - ee) Agr + = (Amq — to] a 8) 
Where 

Aqs qos 
Amq = Xml [= + Sara EE) 

Ads Adr 
Amd = Xl [oe en teeter ee (10) 

Lim sas xan 


Then substituting the values of the flux linkages to find the currents; 


Iqs = = [AGS = DOG | os sncdv ine ehacentoinatonm saan rdaeenad ecole) 
Ids = m5 [Ads — AMA]. .0.....-- sscsnennnennnnn vevuveceecesse (13) 
iis, + [Agr = NAG a pcre tant taeteadtn china healt esnareninnsie tA) 
Idr = = AGH = AMG |. x icvavwieacowntitiimbarnintanamanpaidesanens (1S) 


Torque and Rotor Speed Equation 
Based on the above equations, the torque and rotor speed can be determined as follows [14]: 


3 /P 
Te= ’ (=) — [AdsIqs — AGSIAS]....0. ooo. ooo. esssssssctssssssseessssneer ve es (16) 
wr = fz (re Ne, dae eree ace eet) 


For squirrel cage induction motor, the rotor voltages Vgr and Vdr in the flux equations are set to zero since the rotor cage 


bars are shorted [1]. The symbols as used above can be defined as follows; 


P= number of poles 

a@b = motor angular electrical base frequency, 
or =rotor angular electrical speed, 

Rs =stator resistance, 

XIs = stator reactance, 

Rr = rotor referred resistance, 

Xlr = rotor referred reactance, 


Xm = magnetizing reactance, 
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Xml = machine equivalent stator reactance, 

T,, =electromechanical torque, 

TL = externally applied mechanical load torque, 

J =Moment of inertia, 

Vas = q-axis stator voltage, 

Vds = d-axis stator voltage, 

Vdr = d-axis referred rotor voltage, 

Var = q-axis referred rotor voltage, 

Iqs = q-axis stator current, 

Ids = d-axis stator current, 

Iqr = q-axis referred rotor current, 

Idr = d-axis referred rotor current 

Simulink Implementation of the Induction Motor Sub-Models 
The derived equations represented above are modeled using the subsystem approach in Simulink. However, the sub-models of the 


various subsystems in the complete Simulink model are shown below; 


Voltage transformation Sub-model 


Figure 3 Clark's transformation 


Valpha_in 


Constant Integrator 


Vbeta_in 


Product3 


Figure 4 Parks transformation 
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Flux linkage Sub-model 


Fmd_in Sub_Fdr 


Figure 5 Flux linkage sub-model 


(i) 


Gain4 


(ii) 


Gain2 — Integrator2} 


www.discoveryjournals.org OPEN ACCESS 


baled 84 


ARTICLE 


(iii) 


Gain2 —_Integrator2} 


(iv) 


Figure 6 (i)-(iv) represent implementation of the equations (5)-(8) respectively 


Current sub-model 


lqs - 
" Alpha_beta inverse_abc Ic 


Figure 8 This block implements the inverse park and Clarke transformation. 
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Motor speed and Electro-mechanical torque Sub-model 


Gain1 


- Product1 
Ids_in TL 


Integrator 


Figure 9 The implementation of the torque equation Te (16) and the angular speed equation ur (17) 


Complete Simulink block Implementation of the three phase induction motor 


a 
V 
a Goto Fmq Constant 


Goto Fmd 


4 2 
ro sa 
Current_subsystem Cu 


From wrt 


Flux_Linkage 


TL3 


TL4 


Gain 


Rotor Current 


To Workspace 


«| «OKO 


To Workspace1 


TL 


Cr a! 


Clock 


Figure 10 Complete Simulink block model of three phase induction motor 


Simulation of the model 
The induction motor rated parameters are listed in table 1. 


Table 1 Induction motor rated parameters 


To Workspace? 


Parameter Symbol Value 
Phase Voltage Vp 220V 
Frequency 7 50Hz 

No. of Poles P 4 

Moment of Inertia J 0.089kg.m? 
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Rotor resistance Rr 0.8160 
Stator resistance Rs 0.4350 
Rotor reactance Xlr 0.754Q 
Stator reactance Xls 0.7540 
Magnetizing reactance xm 26.130 
Base angular frequency cob 314.2rad/s 


The induction motor was simulated with variable load stepped torque as shown in table 2. Also, the performance of this 
induction motor with different rotor and stator resistance and reactance values were simulated. Table 3 contains different rotor and 
stator parameters used in this work. 


Table 2 Different motor loading conditions 


Motor Loading TL (N.M) 
No Load 0 

Low Load 100 
Medium Load 550 
Rated Load 1200 


Table 3 Values for motor variable parameters 


Variable Parameters 

Parameters Resistance (Q) Reactance (Q) 
Low stator 0.16 0.016 

High stator 1.0 2.1 

Low rotor 0.1 0.016 

High rotor 2.0 2.1 


A graph of Electromagnetic torque (N.M) Against time (s) 
Rated Parameters 
2000 T T 


1 2 3 4 5 6 


High Stator Resistance 


1000 |- 4 
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= 
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Low Stator Resistance 
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Te(N.M) 


Te(N.M) 


Te(N.M) 


Te(N.M) 


High Rotor Resistance 
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Figure 11 Graph of Electromagnetic torque (N.m) against time (s) 


A graph of Motor Slip Against time (s) 
Motor Rated parameters 
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Low Stator Resistance 
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High Rotor Resistance 
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Figure 12 Graph of Motor Slip against time (s) 


2. DISCUSSION OF SIMULATION RESULTS 


From the above simulation results, overshoot in electromagnetic torque was observed across all parameter variations in the motor 
as shown in figure 11 although it occurred higher in the cases of low stator resistance, high stator and rotor reactance, and high 
rotor resistance parameter values. Precisely, a very high overshoot was observed in motor torque both in starting and load step 
inputs at low stator resistance. This resulted to the motor slip spanning between values of 1 and -2.3 as shown in figure 12. This 
overshoot is about four times its rated value. However, large increase in motor slip due to applied load input was observed with 
either high stator or rotor resistance values. Also, the motor experienced jerky start with high stator resistance value. The motor 
transient lasted for a longer time in low stator resistance than in high stator resistance value. A low overshoot was also observed in 
this high stator resistance value. Transient period at low reactance lasted longer compared to that of the rated and high reactance 
values.That notwithstanding, both stator and rotor reactance increase are capable of causing jerky motor start amidst high 
overshoot at load step inputs. 


3. CONCLUSION 


A 126-hp rated, three-phase induction motor with stepped load at different levels of ONm, 100Nm, 550Nm and 1200Nm was 
modelled in this work by the use of Simulink blocks in Matlab. The motor stator and rotor parameters were varied for motor 
performance studies. From the simulation results, during motor start-up, loading and unloading conditions, the asynchronous motor 
draws large currents, produce voltage dips, oscillatory torque. These results vary with different motor parameters. Low stator 
resistance, high motor reactance, and high rotor resistance parameter values cause high overshoot in motor torque and speed. 
However, low stator resistance causes increased motor transient time during motor start and step inputs. Also, High rotor resistance 
causes wide variations in speed during load step inputs. Increased motor reactance, whether rotor or stator, causes high overshoot, 
jerky start and decreased power factor. Following the observed characteristics it would be very pertinent to recommend that; (i) Soft 
starting of the induction motor should be adopted to avoid voltage dips and machine breakdown. This can be achieved through 
induction motor star-delta starting technique. (ii) Enough vents and cooling system should be provided in the motor to decrease the 
temperature of the motor during operation. Hence, preventing rise in stator resistance. (iii) Due to the wide variations in speed 
during load step inputs, external rotor resistance which is often used to achieve maximum torque at different slip should be 
disconnected during motor operation for efficient and good motor torque speed characteristics. (iv) In variable frequency drive 
(VFD), the operating frequency to rated frequency ratio should be maintained at a value less than unity during motor starting. This 
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operating state will cause a decrease in stator and rotor reactance. Hence, reduced jerky motor start with decreased overshoot and 


improved power factor is achieved. 
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